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Abstract 

1). 

For charged black hole, within the grand canonical ensemble, the decay rate from thermal 
AdS to the black hole at a fixed high temperature increases with the chemical potential. We 
check that this feature is well captured by a phenomenological matrix model expected to describe 
its strongly coupled dual. This comparison is made by explicitly constructing the kink and 
bounce solutions around the de-confinement transition and evaluating the matrix model effective 
potential on the solutions. 
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1 Introduction and Summary 



Quite a while back, Hawking and Page (HP) argued that at temperature below To = y/2/(nl), 
thermal AdS backgrouncQ with a non-contractible Euclidean time circle j3 = 1/T, is stable against 
collapse to form a black holepQ. However, at a temperature higher than Tq, two Schwarzschild 
black holes nucleate, which following the literature, we call small and big black holes. While 
the small black hole is locally unstable, the other one has a positive specific heat. In fact, at 
a temperature T greater than T\ = 3/(nl) > Tq, the large black hole dominates the canonical 
ensemble for quantum gravity defined as a path integral over the metric and all other fields 
asymptotic to AdS with the periodicity of the time direction being /3. It is believed that the 
small unstable black hole facilitates this transition. In the dual gauge theory, this represents a 
de-confining transition for M = 4, large N, SU(N) strongly coupled Yang-Mills on 5" 3 . Though 
for gauge theory at finite coupling, there is no systematic way to compute an effective action, 
one may hope to construct a phenomenological model which captures this transition. Indeed in 
[2] , partition function of YM theory was expressed as a matrix integral over the effective action 
involving Wilson-Polyakov loop. This model, which we review shortly, also predicts certain 
qualitative features that are not visible in the bulk supergravity approximation. Our aim in this 
paper is to put this matrix model to further qualitative checks. This is described in the following 
paragraph. 

As we discussed in the beginning, HP transition is expected to be facilitated by the small 
unstable black hole. A particularly simple way to qualitatively analyze this transition would be 
to construct a bulk effective potential and to find instanton and bounce solutions at and above 
the HP temperature. In this regard, it turns out that the Bragg- Williams (BW) prescription is a 
particularly suitable framework (see [3] for details). In this approach, one constructs an approxi- 
mate expression for the off-shell free energy function in terms of the order parameter and uses the 
condition that its equilibrium value minimizes the free energy. For describing HP transition from 
thermal AdS to the big black hole, horizon radius is an order parameter of the effective potential 
and the temperature then acts as an external tunable parameter. At the saddle point of this po- 
tential, however, the temperature gets related with the horizon radius in a way that is consistent 
with black hole thermodynamics. As we will describe shortly, following Coleman's prescription 
|4J, it is particularly simple to construct an exact instanton at HP temperature T = T\ and the 
bounce above T\ - a region where the big black hole becomes globally stable. Subsequently, one 
can compute the semi-classical decay rate from one phase to the other. Introduction of electrical 
charge into the system makes the chemical potential non-zero. The black hole in question is now 
the Reissner-Nordstrom black hole. In the grand canonical ensemble, as it was pointed out in [5], 
provided the chemical potential is less than a certain critical value, HP transition still persists. 
However, the critical temperature now gets a dependence on the chemical potential. It turns out 
that an appropriate BW potential can be constructed for which both temperature and chemical 
potential acts as independent tuning parameters. Exact instanton and the bounce solutions are 
easily computable. While computing these, we find that decay rate from AdS to the black holes 
phase, perhaps not surprisingly, increases with the chemical potential. Introduction of charge, 
in the boundary, is represented by turning on the R-charge of N = 4 SYM. As we discuss later, 
in the grand canonical ensemble, various qualitative features of this gauge theory with non-zero 

*In our notation I is related to the AdS scale. Later, we work with units such that I = 1. 
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R-charge can still be incorporated into the earlier matrix model provided we make appropriate 
modifications. We find that the instanton and bounce solutions across the de-confining transition 
within this modified model can be semi-analytically calculated. Subsequently, our analysis also 
shows that, as expected from the bulk considerations, there is an enhancement of the decay rate 
with the chemical potential. 

The paper is structured as follows. In the second section, we consider the gauge theory with 
zero chemical potential. In the bulk, we introduce the BW potential and discuss various kink 
and bounce solutions across the two phases and compute the decay rate. After a brief summary 
of the matrix model of [2], we then compute the kink, bounce solutions and decay rate around 
the de-confining transition. Section three addresses similar issues in the presence of a chemical 
potential in the grand canonical ensemble. Here we minimally modify the above matrix model to 
incorporate the effects due to the chemical potential. Once this is done, we go on to check other 
aspects of this model. In particular, we find that, as predicted from the bulk analysis, the decay 
rate from confining to the de-confining phase increases with the chemical potential. 



2 Black hole and their matrix dual: zero chemical potential 



In what follows, we first consider the gravity side, analytically compute instanton and bounce 
at and above the HP temperature. We use this to compute the value of the effective action 
on the solution. This provides us with the decay rate in the bulk. We then briefly review its 
phenomenological matrix dual and carry out similar computations across the de-confinement 
transition within this matrix model. 



2.1 Bulk 



For the Schwarzschild-AdS black holes, the BW potential is given by [3] 

T(r,T) = E-TS = 3(r 2 + l)r 2 -AirTr 3 . (1) 

Here r measures the horizon size of the black hole, while T, E, S are temperature, energy and 
entropy of the hole. Note that all the quantities are made dimensionless by multiplying them 
with appropriate powers of AdS radius /. / does not appear explicitly as it has been set to 1. 
Furthermore, we have set the Newton's constant also to unity. T in (pQ) should be understood as 
a tuning parameter whose dependence on r follows from the on-shell equation 

The thermal AdS potential can be obtained from (Q]) by setting r = 0. Further, the first order 
transition arises when T = and equation ([2]) is satisfied. The solution of these two equations 
are 

3 

r = 1, and T x = — . (3) 

Z7T 
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In figure (Pa), we have plotted ([I]) for various temperatures around the transition point. While 
stable minimum shown by the dotted line represents the stable big black hole for T > Ti, the 
stable minimum at r = 0, shown in the dashed line, represents thermal AdS. 

Following f4], we need to invert the potential, and, look for a kink solution connecting 
degenerate minima at r = and r = 1 at T = T\ shown by the solid line in ([T]b). On the other 
hand, for T > Ti, we need to find the bounce solution connecting meta-stable and stable minima 
of the dotted line in (Q]b). 




(a) (b) 



Figure 1: (a) This is J-(r,T) plotted as a function of r. The dashed, solid and dotted lines are 
for T < Ti, T = T\ and T > T% respectively, (b) Inverse of the potential plotted as a function of 
r. 



This is done by solving the equation 

£ = y/ZF(r,T), (4) 

where r is the Euclidean time. This can be easily integrated with right boundary condition 
leading to the solution 

ggV^r 

r = (5) 

1 + 4irTe^ T + (4ir 2 T 2 - 9)e 2 ^ ' 

This solution is plotted in figure ([2]) for T >T%. For T = T%, it gives a kink solution connecting 
the minima at r = and r = 1 representing the thermal AdS and the black hole phase. At a 
temperature above Ti, we get a bounce corresponding to the semi-classical tunneling from the 
metastable AdS phase to the stable black hole phase. 
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(a) 

Figure 2: Figure shows the instanton and the bounce solutions. The solid line is for T = T% 
3/(2tt) and the dotted line is for T = 3.15/(27r) which is greater than the critical temperature. 



It is straightforward to compute the Euclidean action on ([5]). We get 



S E = J pF{r,T)dr 

-9 + 6vr 2 T 2 + 7iT(-9 + 4vr 2 T 2 ) log (3 - 2vrT) 



V2(- 



9^3 



1 

2T 1 



-3 + 27rr)[3V6-4vrr(2 + vrT) 



+\/3vrT(3 + 2vrT) log (3 - 2vrT + vW3 - 2vrT)] 
For T = T\, the Action evaluates to 



1 



The tunneling rate is therefore, 



_ 1 



(6) 
(7) 
(8) 



2.2 Boundary 



Having analyzed semi-classical decay of thermal AdS to a stable black hole at high temperature, 
we now turn our attention to the boundary dual. As mentioned previously, the framework within 
which we will search for similar kink and bounce solutions is the one proposed in pj^l . It is a 
phenomenological model where partition function of YM theory was expressed as a matrix integral 
over the effective action involving Wilson-Polyakov loop (trU)/N. It contains two parameters 
a and b which, at fixed 't Hooft coupling, depend on the temperature. The partition function is 
given by 



Z(T) 



dU 



a{T) till tr[/t + trU trl7+) 5 



(9) 



^This model was further studied and extended in [6] [51 [9] . 
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Figure 3: Plots of a(T,0) and b(T,0). This is taken from [7J 



This model was studied in some detail in the original paper [2] and was further reviewed in section 
5 of [7J. We will be very brief here in describing the model. The effective potential arising from 
this model can be written in terms of p = trU trW /N 2 as 

V( P ) =^p 2 - b - P * forO<p<i 

= ~\ P 2 ~\p'-\ log [2(1 -P)} + \ f <* \ < P < 1. (10) 
The saddle point equation is given by 

* 1 
ap + bp i = p for > p > - 

= — !— for - > p > 1. (11) 
1 - p 2 ~ p ~ v 7 

The restrictions on a and 6 are a < 1,6 > 0. We note that the comparison between this model 
and the bulk action is valid only when one neglects the string loop effects. The temperature at 
which the supergravity description breaks down is identified with the Gross- Witten transition in 
the matrix model. 

The temperature dependence of the parameters can be fixed as follows. Let us consider a 
temperature T > To, To being the black hole pair-nucleation temperature. At this temperature, 
we can writ^l 

2ap\ 2 + 2bp\ 2 + log (1 - pi, 2 ) + / = -J 1)2 , (12) 

where are the actions for the large and the small black holes. p\ i are the corresponding 
solutions in the matrix model. The constant / = log (2) — 1/2 is introduced to make the potential 
continuous at p = 1/2. We further note that pi j2 also have to satisfy the saddle point equations 
(fTT|) . These four equations uniquely determine a,b and px >2 at a specific temperature. Further, 
continuing this process for various other temperatures, we get the dependence of temperature on 
a(T) and b(T). This can be done numerically and the result is shown is shown in figure ([3]). 

We now proceed to look at instanton and bounce solution of this model. Though the 
method is same as the previous one, in order to proceed, we face a problem. The minimum at 

^For all the numerical calculations, we set I, the Newton's constant, volume of the three sphere to 1. 
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non-zero p always appears in a region p > 1/2. At the de-confining temperature T = T\ = 3/(2?r) 
this minimum appears at p = p\ = 0.856, a = 0.316,6 = 1.193 and degenerates with the one at 
p = dl • So to get the complete instanton solution, we need to include integration of V in (jlOD in 
the second region. However, because of the presence of the log term, we are unable to do an exact 
integration. In this region, we therefore resort to numerical method. In the region < p < 1/2, 
the kink solution is given by 

p(r)= 4 ( 1 ~ a ) e l ~ a J_ for0<p<-, (13) 
PV ' 4 + 6(1 - a) tfVttr ~ H - 2' v ' 

for a = 0.316,6 = 1.193. For the rest of the region, we numerically integrate the potential. 
The result is shown in figure dU). For T = T\, a = 0.316 and 6 = 1.193 the critical point is at 
P = Po = 0.856. Similar to the bulk case a calculation of the action for the instanton solution 
p will give us the width of tunneling between the two solutions. A numerical integration of the 
instanton action evaluates to 

rPn i 

2V(p)dp = 0.175 (14) 



Hence the tunneling rate becomes 

0.175 

r = e — — . (15) 

For T > T\ the p = solution becomes metastable and, similar to the bulk case, one obtain 
bounce solutions. This is shown in figure ©. 



^ These values follow from [7]; see figure 8 there. 
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0.9-1 




Figure 5: (a) This is V(p) plotted as a function of p for T > Tj, a = 0.316 and b = 1.220. (b) 
The bounce solution for the same potential with the turning point po = 0.810 



3 Black hole and their matrix dual: non-zero chemical potential 



In this section, we will continue similar exploration in the presence of non-zero chemical potential. 
Our whole analysis will be carried in the grand canonical ensemble. We first study this problem 
in the bulk. In the boundary, we generalize the matrix model to include the effect of chemical 
potential. We then find the kink and bounces in this model. Calculating the decay rate from 
confining to the de-confining phase, we show that rate increases with the chemical potential. This 
is consistent with the expectation that arise from the bulk behavior. 



3.1 Bulk 



In the presence of electrical charge, we need to consider the Reissner-Nordstrom black holes. The 
BW free energy is easily computed and is given by [3] 

F(r, T, p) = E-TS-Qp 

= 3r 2 (l - -p 2 ) -4vrr 3 T + 3r 4 . (16) 
3 

where Q is the charge and p is its conjugate chemical potential. We will restrict ourselves to the 
case where p < V3/2. At the saddle point, 



3vrT + V9vr 2 T 2 - 3(6 - 8p 2 ) 
r= 7. • (17) 
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This can be inverted to get black hole temperature. While reducing the temperature below a 
critical value (Ti), this hole becomes globally unstable and system crosses over to the thermal 
AdS phase. This critical temperature depends upon the chemical potential and is given by 



3 Jl - Im 2 



Tl = (18) 



As in the previous section, the kink solution within this effective description can be con- 
structed at the critical temperature. It is given by 



r(Y) 



2(3 - 4/i 2 ) e v / 3 Z ^ 



1 + 2^9 - 12^ 2 eV 3 -V * 



(19) 



That the solution indeed interpolates between two phases, namely, the thermal AdS and the black 
hole, can be seen from the figure ([6h.). For T > T\, the thermal AdS phase becomes metastable 
and we expect to find a bounce solution connecting this metastable phase to the stable one 
representing the black hole. Indeed this bounce can be explicitly constructed and is given by 



r t 



2(3 - V)e 



1 + 4eV3-4M 2 r nT + e 2V3"4M 2 r ( 12fJ 2 + 4 ^2 T 2 _ g) 

This is plotted for different values of \x in (Jib). 



(20) 





(a) 



(b) 



Figure 6: (a) Instanton solution at T = T\ for \i = .5, .3 and .1 are shown by solid, dashed and 
dotted lines respectively (b) Bounce solution for T > T\ for different \i. Solid, dashed and dotted 
lines are for fi = .5, .3 and .1 at temperatures .4, .45 and .48 respectively. T\ for these \i values 
are 0.390, 0.447, 0.474 respectively. 



It is straightforward to compute the Euclidean action on the instanton solution (|19p. We 



get 



Se 



4/<- 



2T(r,T,n) 
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^2(^(3^3- V(_ 3 + 2vr 2 T 2 + 4/?) 
+\/3vrT(-9 + 4vr 2 T 2 12// 2 )Loc/[-27rr + a/9 - 12^ 2 ]) 



+(>/2(-3 + 4^ 2 )(-12tt 3 T 3 + 47r 2 T 2 ^/9 - 12^ 2 + 15ttT(3 - 4/i 2 ) 

+6^9 - 12/i 2 (-3 + 4/i 2 )) + ttTJ9- 12/x 2 (9 - 4vr 2 T 2 - 12^ 2 ) J 9 - 12^ 2 - 2tttJ9- 12 fi 2 



log(-2vrT + a/9 - 12^ 2 + \/ 18 - 24,u 2 - 47TTJ9 - 12/i 2 )) 



/(27^ (-3 + 4/i 2 ) + 12/i 2 + 2irTy/9 - 12^ 2 ) )) 
At T = T\ given in (fTSI) . we get 



(3-4/. 2 
9^/2 



The tunneling rate is therefore, 

We see that decay rate increases with \x. 



(21) 

(22) 
(23) 



3.2 Boundary 

We now generalize the matrix model considered in the previous section to include the effect of the 
chemical potential. It turns out that minimal modification requires one to keep the structure of 
the potential same as Q, but to make a and b dependent on temperature as well as the chemical 
potential. One can now follow the same strategy as before to figure out the dependence of a, b 
on n and T. The only difference is that now in (|12p . we need to substitute Ii2 for the Reissner- 
Nordstrom black hole. Carrying out the computation numerically we get a, b as shown in the 
figure ([7]). A clear pattern that we see from here is that both these two parameters increase as 
we increase the chemical potential. 



0.36 0.38 



0.42 0.44 
Temperature (T) 



0.46 0.43 



0.42 0.44 
Temperature (T) 



0.46 0.43 



(a) 



(b) 



Figure 7: Plots show the behaviours of a(/x, T) and 6(/i, T). 
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Once we have the data, we can go ahead and numerically compute the instanton, bounce 
and the decay rate. This is similar to the previous section with a little more complicacy due to 
non-zero \x. We skip the details and show the results of our computation in the following figures. 

0.9 -i 0.8 -i 




-i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

-4-3-2-10 1 2 3 4 




I ' 1 1 1 1 1 1 1 1 1 1 1 

-6 -4 -2 2 4 6 



% X 

(a) (b) 

Figure 8: Boundary instanton and bounce curves. Compare with the previous figure 
dH) for the bulk case, (a) Solid, dashed and dotted curves correspond to instanton solutions 
for Ou,Ti) = (0.5,0.390) (0.3,0.447) (0.1,0.474) respectively, (b) Bounce solutions at T > 
T\ for different potentials. The solid, dashed and dotted curves correspond to (/i, T) = 
(0.5,0.4), (0.3,0.454), (0.1,0.48) respectively. 



The effective action on the instanton solution can be evaluated numerically for this matrix 
model at different chemical potentials. The result is shown in figure (9). It is satisfying to 
find that, as in the bluk computation, the value of boundary effective action decreases with the 
increase of the chemical potential. This leads to an increase in the tunneling rate. We also 
note that the boundray dual is a phenomenonological model and is expected to capture only the 
qualitative features of the gravity dual. The functional forms of J-{r) and V(p) being different, 
it is not surprising that the two lines in figure (9) do not coincide. 

In our study, we have adopted an effective approach towards the bulk gravity. Though this 
approach is particularly simple, enabling us to analytically construct kinks and bounce solutions, 
it would be instructive to directly analyze the gravitational bounce, compute the action against 
the bounce and then compare the results with the matrix model results. We leave this for our 
future study. 
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Figure 9: Behavior of effective matrix model potential calculated on the instanton is shown by 
the dotted line for different chemical potential. For comparison, the corresponding bulk action is 
also shown in the plot (dashed line). In both cases, actions decrease with fi. 



Acknowledgments: B.C. would like to thank Institute of Physics, Bhubaneswar for hospitality 
during the course of this work. 



References 

[1] S. W. Hawking and D. N. Page, Commun. Math. Phys. 87, 577 (1983). 

[2] L. Alvarez-Gaume, C. Gomez, H. Liu, S. Wadia, Phys. Rev. D71, 124023 (2005). [hep- 
th/0502227]. 

[3] S. Banerjee, S. K. Chakrabarti, S. Mukherji, B. Panda, Int. J. Mod. Phys. A26, 3469-3489 
(2011). [arXiv:1012.3256 [hep-th]]. 

[4] S. R. Coleman, Phys. Rev. D 15, 2929 (1977) [Erratum-ibid. D 16, 1248 (1977)]. 

[5] A. Chamblin, R. Emparan, C. V. Johnson and R. C. Myers, Phys. Rev. D 60, 064018 (1999) 
[arXiv:hep-th/9902170]. 

[6] P. Basu and S. R. Wadia, Phys. Rev. D 73, 045022 (2006) [arXiv:hep-th/0506203]. 

[7] T. K. Dey, S. Mukherji, S. Mukhopadhyay, S. Sarkar, JHEP 0704, 014 (2007). [hep- 
th/0609038]. 



12 



[8] T. K. Dey, S. Mukherji, S. Mukhopadhyay, S. Sarkar, JHEP 0709, 026 (2007). 
[arXiv:0706.3996 [hep-th]]. 

[9] T. K. Dey, S. Mukherji, S. Mukhopadhyay and S. Sarkar, Int. J. Mod. Phys. A 24, 5235 
(2009) [arXiv:0806.4562 [hep-th]]. 



13 



